The influence of spin-flip scattering on the preparation and detection of a single spin 
state in a quantum dot attached to a spin battery 
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Recently, the possibility of an all electrical scheme of preparation and readout for a single spin 
state in a single quantum dot attached to spin biased leads has been shown [F. Chi et al Phys. Rev. 
B 81, 075310 (2010)]. However, spin scattering mechanisms have been omitted. To remedy this 
lack we consider the influence of the spin-flip scattering process on the proposed preparation and 
readout scheme. 
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I. INTRODUCTION 

Coherent manipulation of a single spin state in quan- 
tum dots (QDs) is crucial for spintronics and quantum 
information physics. A qubit, which is the basic require- 
ment for quantum information processing, can be realized 
as an electron's spin on a single-electron quantum dot^. 
Quantum dots seem to be promising elements for spin- 
based quantum computer as their properties are rather 
easily and electrically controllable. Preparation and ma- 
nipulation of a single spin state in a QD has been real- 
ized by various optical and magnetic means^Ti^ or spin to 
charge conversion^' ^ techniques. However, to our knowl- 
edge, there is no experiment allowing all-electrical control 
of a single spin state in a QD (altough such an attempt 
should reduce the complexity of the experiment). 

Recently, Chi et al have proposed a theoretical scheme 
for the preparation and readout of a single spin state 
in a QD utilizing spin battery^ realized in recent 
experiments^!^. Manipulation and detection of a spin 
state have been also investigated in a double quantum 
dot system coupled to the spinii or chargei^iH biased 
leads. However, spin scattering mechanisms in a sin- 
gle QD have not been taken into account^. Intradot 
spin-flip processes can dramatically change a dot's trans- 
port characteristics and lead to suppression of the tun- 
nel magnetoresistance-i^^ii^. At low temperature the two 
most important mechanisms in the QD responsible for 
spin-flip phenomena are spin-orbit interactions^^ (SOI) 
and the hyperfine interactionii (HI) which couples the 
electrons spin to an effective magnetic field induced by 
nuclear spins. However, for B ^ Bnuciear ~ 3mT the 
HI mechanism is suppressedi^ and the spin relaxation 
is then mainly due to spin-orbit interaction. Many low 
temperature experiments have been realized for measur- 
ing the spin relaxation time in semiconductor quantum 
dots§.iiilii^"^. Although all the ingredients of the QD's 
spin-based quantum computer are available, it is still a 
great challenge to build such a device, accessible for com- 
mercial applications. 

The electron spin state can relax due to thermal de- 
cay, which prevents spin-based quantum computation at 
room temperature. Another dot's spin-flip mechanisms 



can be induced by higher order tunneling events. Namely, 
when the dot is in the Coulomb blockade regime, co- 
tunneling processes may flip the spin on a dot^-?. How- 
ever, incoherent processes, like inelastic transitions and 
co-tunneling, provide rather a minor contribution to the 
spin-flip mechanism^. It is also worth mentioning that in 
the case of strong coupling to the leads at low tempera- 
ture significant impact on the spin scattering mechanism 
can originate from Kondo resonance. 

In this paper we analyze the influence of the intradot 
spin-flip mechanism on the control of a single spin state 
in a quantum dot attached to a spin battery. A spin 
battery2ii2iiS^ allows to inject pure spin current into the 

d0t?4. 

The paper is organized as follows. In section |ll] we 
describe the model and theoretical formalism. Numeri- 
cal results are presented and discussed in section Hill A 
summary and final conclusions are gathered in section 
llYl 



II. MODEL AND THEORETICAL FORMALISM 

We consider a single-level quantum dot coupled to one 
spin battery and one normal lead. Then the Hamiltonian 
of the system is of the form: 

= ^ £kQ<TCk„o.CkQo- + ^ f^aqtqa + Un„ng 

kacT CF 

+R{q\q, + qlq,) + Y^T.^^^lo..'!'^ + h'^')- (1) 
ka ly 

The first term describes the left [a = L) and right 
{a — R) lead in the non-interacting quasi-particle ap- 
proximation. Here, cj^Q.^. (ckacr) is the creation (anni- 
hilation) operator of an electron with the wavevector k 
and spin a in the lead a, whereas Ekao- denotes the cor- 
responding single-particle energy. The next three terms 
in the Hamiltonian ([1]) describe the quantum dot. Here, 
n„ — q^qa is the particle number operator (tr =t, 4-), Co- 
is the discrete energy level of the QD and U is the intra- 
dot Coulomb integral. The parameter R in the Hamilto- 
nian (HI) describes the spin-flip transition amplitude. The 
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spin-flip term is assumed to be coherent, in the sense that 
the spin- flip strength R involves reversible transitions be- 
tween up- and down-spin states on the dot. Such an effect 
may originate from a spin-orbit coupling-- in the dot or 
from the transverse component of a local magnetic field 
applied, for example, within the electron spin resonance 
techniques^. The last term of Hamiltonian ([IJ describes 
electron tunneling between the leads and dot, where V" 
are the relevant tunneling matrix elements. Coupling of 
the dots to external leads can be parameterized in terms 
of r"(e) = 2TT\V^\'^Pa, where Pa is the density of states 
in the ath lead. In the wide band approximation F" is 
constant within the electron band, r"(e) = F" = const 



Poo = 

Paa = 

+ 



Pas = 

+ 



P22 = 



In above equations, = + F^ with F^{x) ~ 
J2a^»crf{x-iJ.acr) and F" (x) = I]„FQ^[l-/(a;-^QCT)], 
for X = {eCT±-R, f-a+U ±i?}. Here, f{e—pa<y) is the Fermi- 
Dirac distribution function in the ath lead. The ther- 
mally averaged diagonal elements of the density matrix 
(poo, Pcrcr, P22) refer to the probability of the dot's state 
being empty, singly occupied by an electron with spin a 
and doubly occupied, respectively. In turn, the nondiag- 
onal terms p^^ describe coherent transitions between a 
and a spin states on the dot. Rate equations 1^ together 
with the completeness relation poo + J2a P<^<^ + /^22 = 1 
give us full information about the dynamics of the sys- 



for e G (-VF/2, W/2), and F;^(e) = otherwise. Here, W 
denotes the electron bandwidth. 

To find equations governing the dynamics of the sys- 
tem in the weak coupling regime we follow Ref.^ except 
as regards the derivation of the dot's Green's functions. 
In contrast to Ref.— we do not restrict our model to weak 
spin-fiip strengths by calculating the dot's Green's func- 
tions assuming vanishing R. Here, the dot's Green's func- 
tions are obtained for finite R. Details of the method are 
displayed in the Appendix. For the sake of simplicity we 
assume spin degenerate dot energy level (e^ — = eo). 
The rate equations for dot's expectations values of the 
density matrix elements acquire the following form: 



-R)]ps.} 

- R)] Paa 



(2) 



tem. The spin-dependent dot's occupation numbers are 
expressed in the following way Ua = Pan + P22 ■ 

Current flowing from a lead to the dot is obtained from 
the standard definition: 

Ji = -e{Nc.)^~i\{[H,N^]), (3) 

where is a occupation number operator in a lead. Af- 
ter performing onerous calculations the current formula, 
in the weak coupling approximation, acquires the follow- 
ing form: 



I 

^ Y,{[F-{eo + R) + F-(eo - R)]p.. - [F+{eo + R) + F+(eo - i?)]poo + [P- {e^ + R) - F' {e^ - 

a 

i{[F+(eo + R) + F+(eo - R)\poo - {K i^o + i?) + P- {t^ - R) + F+{eo + U + R) + F+(eo + U 
[F^ {eo + U + R) + F^ {eo + U- i?,)]p22 + 2iRp,^ 

[F-(eo + R)- F-{eo - i?) - F+{eo + U + R) + F+ {e^ + U - R) + 2iR]p,,}, 
i{[F+(eo + R) - F+{eo - i?)]poo - [F-{eo + R) ~ J^"(eo - R) + 4.tR]p^^ 
[F+{eo + U + R)-F+{eo + U + R) + AiR]p„„ 

[F-(eo + R)+ F-{€o -R) + F+ieo + U + R) + F+ {e^ + U - R)]pa^}. 

\ I]{[^.+ (eo + U + R) + F+{eo + U - R)]psa - [F+(eo + U + R) ^ F+ {eo + U - R)]p^s 

a 

\F- {eo + U + R) + F- {eo + U- R}]p22}- 

I 



I 

= ^Y.'^i^^^'^^'o + R) + F+,{eo - R)]poo ~ [F-Aeo + R) + F-^{eo - R)]p,, 

a 

+ [F+^eo + U + R) + F+^eo + U- R)]p^^ - [F-^eo + U + R) + F'^eo + U- R)]p22 

- [F-^eo + R)- F-^eo - R) + F'^eo + U + R) - (eo + U- R)]p^,}, (4) 



where ?R.[A] denotes the real part of A. Finally, current passing through the dot can be symmetrized in the fol- 
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FIG. 1: Dot's occupation numbers as a function of the spin 
bias voltage calculated (a) in the absence of the spin-flip relax- 
ation {R = 0) and for the indicated values of the temperature, 
(b) for the indicated values of parameter R and for tempera- 
ture fesT = 0.01. 



lowing way: J = {Jl — Jr) /2. The spin current is defined 
in the following way: Js — {h/2e){Jf — Jj,), where is 
the spin-cr electron contribution to the charge current. 



III. RESULTS 

In numerical calculations we take the parameters from 
RefH. As the charging energy is the largest energy in the 
QD we choose it as the energy unit, i.e., U = 1. Then, 
other parameters of the system are expressed in units 
of the Coulomb interaction parameter U. We set the 
dot-lead coupling to be = = 0.02. Moreover, we 
assume that only the left lead provides spin bias, whereas 
the right electrode is the normal lead. Specifically, we set 
^]^-^ = eVs, ULi = -eVs and hr^ = ^iri = = 0, where 
Vs is the applied spin bias. No (charge) bias voltage is 
applied. 

Following RefjS we start from analyzing the prepara- 
tion stage in the stationary limit. In this case we set the 
left hand sides of Eqs.(I2|) equal to zero, i.e., Pij{t) — 0. 
Steady state master equations together with the com- 
pleteness relation allow us to obtain relevant expectation 
values of the density matrix elements and dot occupa- 
tion numbers of each spin component (n^ for a ^tii)- 
At the beginning we calculated in the absence of spin- 
flip processes, which is equivalent to making the assump- 
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FIG. 2: Time evolution of the dot's occupation numbers cal- 
culated for initial state |0, 1) and for the indicated value of 
parameter R, and for eo = —0.5, fcsT — 0.01. 



tion that R = 0. Figure [TJa) shows dot's occupation 
numbers as a function of the applied spin bias calcu- 
lated for the dot level situated symmetrically between 
^R and (eo = —eVs/2). In such a prepared sys- 
tem QD tends to be occupied by an electron with spin 
up when a sufficiently large spin bias is applied. With 
increasing temperature one needs to apply greater spin 
bias to obtain a QD in the state with a spin-up electron 
with probability close to unity. This indicates the possi- 
bility of preparation of a single spin state in a quantum 
dot by using a spin bias which has been shown in Refj^. 
Moreover, a sufficiently high spin bias is now available in 
experiments^ii^. However, spin-flip processes may have 
a significant impact on the above picture. Thus, now we 
consider the influence of these processes on the prepara- 
tion stage. In FigHfb) we plot dot occupation numbers 
for different values of the parameter R (proportional to 
the strength of the spin- flip processes) . Generally, when 
R is nonzero the occupation number n-^ (ni) is reduced 
(enhanced). However for small R and sufficiently large 
spin bias {eVs > 0.68) this suppression of n-^ (enhance- 
ment of ni) can be partially recovered. Thus, if the 
intradot spin relaxation is weak enough we are able to 
prepare dot in a given state. On the other hand, when 
the spin-flip processes are strong enough {R > F) the 
dot can be occupied with an electron with spin up or 
down with almost equal probability. Then, even increase 
in the spin bias voltage is not able to enforce the dot's 
occupation by a spin-up electron. In contrast to the 
weak spin-flip case, now, further increase of the spin bias 
leads to increase of both n-f and ni. Moreover, addi- 
tional features in the dot's occupation number depen- 
dence can be noticed. These features are due to splitting 
of the dot's energy level induced by a sufficiently strong 
spin- flip strength i.e., for R ^ T. After the initializa- 
tion stage we lower the dot's energy level in such way 
as to prohibit sequential tunneling events. The sys- 
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FIG. 3: Time dependence of the occupation numbers for different initial states (n|(0), n^(0)) calculated for the indicated values 
of parameter R. Other parameters: to = —0.9, ksT = 0.01. 




FIG. 4: Time evolution of the charge current (Jr) for different initial states (n|(0), 7i^(0)) calculated for the indicated values 
of parameter R. Other parameters: eo — —0.9, fcsT — 0.01. 



tern is then in the Coulomb blockade regime, because 
£0 < ML(T,/^fl and eQ + U > ulut^-b- In order to check 
how long a spin state can be maintained before a flip 
to another state we have derived numerically the rate 
equations ^ and calculated the time dependence of the 
occupation numbers for given initial conditions. In 
the Coulomb blockade the dot is occupied by one elec- 
tron. Thus, poo(O) = = P22(0), and we can write the 
initial dot's state as |n^(0), 714.(0)) = {|1, 0), |0, 1)} or as a 
superposition of spin-up and spin-down states. When the 
intradot spin- flip processes are neglected the state |1,0) 
is very stable, whereas the spin-down state |0, 1) evolves 



into the state |0, 1) after some time^. For a typical quan- 
tum dot at very low temperature this time can be long 
enough for performing many qubit operations. However, 
the spin-flip effect may greatly shorten this time scale. 
In Figl5]we plot the time evolution of the dot's occupa- 
tion numbers for nonzero R and for the dot's initial state 
0, 1). One can notice that even weak spin- flip processes 
have a significant impact on these quantities and hence 
on the spin lifetime. We have also estimated the spin 
lifetime for different values of the parameter R assuming 
U = IQmeV, ksT = O.OIU (T - IK). The results ob- 
tained are shown in TablH When the spin-flip strength is 
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TABLE I: Calculated spin lifetime corresponding to the given 
value of parameter R 
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larger than R > 0.01 (O.lmeF for the assumed parame- 
ters) the spin hfetime is too short for qubit manipulation 
even using (picosecond) optical methods. The time evo- 
lution of each component of the occupation numbers is 
a periodic function of the time. At the times T = 
{n = 1,2,...) the initial dot's state is restored. When 
the dot is initially occupied by a spin-up electron the 
time evolution is analogous to that when the initial state 
is |0, 1) after making the transition n-f-(0, 1) — >■ n^(l,0) 
and 714,(0,1) — >■ n^(l,0). With increasing R the time 
evolution of the occupation numbers is qualitatively the 
same, however, the time scale is correspondingly changed 
(shorter) . 

Finally, we examine the influence of the spin-flip scat- 
tering on the read out process. After performing manip- 
ulations on a qubit state the dot's energy level is further 
lowered until /i^^ > eo -|- J7 > /i/j. Now, spin- up elec- 
trons can take part in the transport. Regardless of the 
spin relaxation processes, when the dot's level is initially 
occupied by a spin-down electron, an electron with spin 
up can tunnel through the dot (onto the dot from the 
left lead and further to the right lead) and a finite cur- 
rent arises. However, when a spin-up electron resides in 
the QD, the spin-down electron can not tunnel onto the 
dot because + U > ^iLi^ whereas the spin- up electron 
can not leave the dot because eo is beneath the transport 
window. As a result the current is suppressed. That is 
why measurements of the charge current should give us 
information about the dot's state. 

When i? = a finite current can be observed when dot 
was initially occupied by a spin-down electron, whereas 
the current vanishes for the state |1,0)>^. After a suffi- 
ciently long time the state |0,1) finally evolves into the 
state |1, 0) and the charge current is suppressed. However 
the above scenario ceases to hold when spin-flip scatter- 
ing is included. Due to the spin-flip process, now, the ini- 
tial dot's state |1, 0) (|0, 1)) can be switched to the other 
spin state and the current can be partially enhanced (sup- 
pressed) as is shown in Fig]?] This switching mechanism 
leads to oscillations in occupation numbers (see Figl3]) as 
well as in the current with period inversely proportional 
to the strength of the spin-flip processes. These oscilla- 
tions are dumped due to the dot-lead coupling and the 
current saturates to a certain value. 



IV. CONCLUSIONS 

In conclusion, we have studied coherent transport 
through a single quantum dot subjected to spin bias in 
the presence of intradot spin-flip processes. We have 



found that sufficiently large spin-flip strength can pre- 
vent preparation and readout dot's spin state by means 
of the present techniques. However, recent experiments 
have shown that the spin lifetime can reach millisecond 
or second time scales^^ii^ which is enough for perform- 
ing many spin operations. That is why, along with the 
authors of Ref.lsl, we believe that a single spin state in 
the (all electrically controllable) quantum dot is a good 
candidate for being a qubit. 
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Appendix: Rate Equations 

To obtain the rate equations ^ we adopt the formal- 
ism presented in Ref.^*". Specifically, we express the dot's 
operator in terms of Hubbard operators: 

q,^\Q){a\+a\a){2\ (5) 

represented by four possible electron states in each dot 
which satisfy the following completeness relations: 

|0)(0|+^H(a| + |2)(2| = i. (6) 

a 

Furthermore, the set of auxiliary operators is introduced 
and the dot operators are expressed by means of the 
slave-boson and pseudofermion operators'^^. 

From the definitions of the Dirac brackets one is able 
to find the commutation (and anticommutation) rules for 
new operators 1^ In the slave-boson representation the 
Hamiltonian of the system ^ acquires the form 

H = ^ekaacL<7Cka<7 + ^eff(/]/<T + d^fi) 

+ R{flf^ + Il.c.) + UdU 

+ EE[C4..(eV. + ^/^d) + H.c.]. (7) 

ka cr 

Here, is the slave-boson operator which creates an 
empty state in the dot, /t is a peudofermion operator 
which creates a singly occupied state with an electron 
with spin a, whereas creates a doubly occupied state 
with an electron with spin a and another electron with 
spin a in the dot. In the slave-particle representation the 
density matrix elements are written in the following way: 
Poo — e^e, Pctct = /J/cr, P22 — d)d- Here, the statistical 
expectations of the density matrix elements {pnn = {Pnn) 
with n = 0, cr, 2) give the occupation probabilities of the 
quantum dot being empty, singly occupied by an electron 
with spin cr, and doubly occupied, respectively. 
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To derive the rate equations we start from the von 
Neumann equation for the density matrix operator; 

'p = i[H,pl (8) 

where p = {poo, Paa, Pa-s-, P22)^ ■ Furthermore, the aver- 
aged equations for density matrix elements are expressed 
by means of dot-lead Green functions, and using the Lan- 
greth theorem^^, the dot-lead Green functions can be set 
down by means of dot's Green functions and free leads' 



Green function s^ ' . The Green functions of the dot, in 
time space, are defined in the following way G^a-'(t,t') = 
{{q.mlit'))) = {{eHt)Ut)\fl,it')eit'))) + 

aa'{{fl{t)d{t)\S{t')f,.{t'))) = Geaa'it,t') + Gd,,'{t,t'). 

Other parts of Gaa' {t, t') vanish for t' = t, and thus are 
omitted as we are interested in the t' — t case. Now, we 
are able to express the rate equations by means of the 
dot's Green's functions and these can be written in the 
following form: 



POO = /'dc.^[crHG>,M + r^(i-rH)G<_H] 

-I- iR{pcrS - Pscr) 

p.A - ^ /dc.5]{(c + r^)[rHG>,H + (i-rH)G<,H-rHG,^^^ 

OL 

+ iR{Paa - Paa) 

= hi ^^E[r"r(^)G^..(^)+r^(i-rM)G^.*M] (9) 

I 

Now, we only need to know the Green's function of the from the corresponding equation of motion for the dot's 
dot. This is derived in the weak coupling approximation operators. Thus, we obtained 





= iTTpaa{S[uJ - 


(eo + 


R)] +S[oj 




-R)]} + 


inpaa{S[u} - (eo + R)] 




= ~~iTrpoa{S[uj 




+ R)] + S[lj - (eo 


-R)]}, 




G%-M 


= iTTp22{S[uJ - 


(eo + 


U + R)] + 


S[iu - 


(eo + f/ 


-R)]}, 




= -iTTpff^{d[LLl 


-(eo 


+ U + R)] 


+ S[u} 


- (eo + 


U - R)]} + i7rpaa{S[uj 


Gf.-A^) 


= iT:paa{S[ljJ - 


(eo + 


R)] +S[oj 


-{eo- 


-^)]} + 


iTTpaa{S[l^ - (co -|- R)] 




= -TTpaoi^i'^ - 


-(eoH 


^ R)] - S[lu 


-{eo 


~R)]}, 




G^,A^) 


= i'!Tp22{S[u} - 


(eo + 


U + R)]- 


S[oj- 




~R)]}, 


G^M 


= ITT Paa {6[UJ - 


(eo + 


U + R)] + 


5[uj- 




- R)]} - iTrpaa{S[uj - 



(10) 



To derive these Green functions we assumed no coupling 
to the leads {V" = 0) and that the leads are taken to be 
in local thermal equilibrium. In contrast to RefsJ^i^li^ 
we did not assume vanishing R during the calculation 
of the above Green's functions. Thus, our equations are 



not limited to small values of the spin-flip strengths (i? <C 
r). Moreover, the form of these Green's functions clearly 
shows that spin- flip processes lead to splitting of the dot's 
energy level. Finally, connecting Eqs. ([T(Il) with Eqs.® 
we arrive at the coupled set of differential equations ^ . 
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